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1. Introduction 

The first law of motion to be analyzed is the Sedov-Taylor solution 

where E is the energy injected into the process, p is the density of the surrounding medium, and t is the time, 
see [1-3]. This equation allows of deducing the energy of the first explosion of the atomic bomb, the Trinity test in 
New Mexico in 1945, which is E fa 10 21 ergs. A second application of this equation belongs to astrophysics and is 
the supernova explosion (SN) in which E ~ 10 ergs, see [4]. A different approach as given by [5] and [6] analyzes 

two self-similar solutions with varying inverse power law exponents for the density profile of the advancing matter, 

_ _ "~ 3 

R , and ambient medium, R s . The previous assumptions give a law of motion R oc t ™- s when n > 5. Another 

example is an analytical solution suggested by [7] where the radius-time relationship is regulated by the decrease 

in density: as an example, a density proportional to R~ 9 gives R oc t 2 ^ 3 : this mechanism as been applied to 

young supernova remnants, see [8]. The case of spherically symmetric thin shells has been recently analyzed both 
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in special relativity, see [9] with application to the mass inflation phenomenon, and in general relativity, see [10] 
where the solutions are split into "ordinary" matter and "exotic" matter. The actual theoretical situation leaves 
a series of questions unanswered or only partially answered: 

• Is it possible to deduce an equation of motion for an expanding shell by assuming that only a fraction of 
the mass enclosed in the advancing sphere is absorbed in the thin layer? 

• Can this new adjustable law of motion for an expanding shell be calibrated on the observed data? 

• Can we deduce an asymptotic behavior for the advancing radius of the thin layer of the type oc t a , where t 
is the time and a an observational exponent? 

• Can we model the momentum carried away by the photons? 

• Can we derive a relativistic law of motion for the expanding shell under the previously outlined hypothesis? 

In order to answer these questions, Section 2 reviews the existing situation and reports a new classical law of 
motion which introduces the two concepts of porosity and photon's losses. Section 3 reviews a first relativistic 
law of motion and introduces a new relativistic law of motion which includes the porosity. 

2. Classical thin layer approximation 

This Section reviews the general purpose power law model, the standard thin layer approximation with constant 
density and the thin layer approximation in a medium with density regulated by a power law. A first new equation 
of motion is introduced which works in a medium with constant density but swept mass regulated by a parameter 
called porosity. A second new equation of motion includes porosity and photon's losses. 

2.1. The power law solution 

The equation of the expansion of a shell can be modeled by a power law of the type 

to 

where R is the radius of the expansion, t is the time, Ro is the radius at t — to, and a is an exponent which can 
be found from a numerical analysis. 
The velocity is 

v(t) = aRoi^Tt^ . (3) 
to 

As an example, Fig. 1 reports the fit of SN 1993J where the radius is growing more slowly than the free expansion, 
R ix t, and more quickly than the Sedov-Taylor solution, R oc t 0A , see equation (1). 
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Radius of SN1993J with power law fit 




Figure 1. Theoretical radius as given by the power law fit represented by equation (2) with a —0.828, Rq— 0.0087 pc and £q — 
0.498 yr. The astronomical data of SN 1993J arc represented by vertical error bars and are extracted from Table 1 in 

[ii]. 

2.2. Standard case 

The thin layer approximation assumes that all the swept-up gas accumulates infinitely in a thin shell just after 
the shock front. The conservation of the radial momentum requires that 



4 3 . 4 3 . 
-nR pR = pRo , 



(4) 



where R and R are the radius and the velocity of the advancing shock, p the density of the ambient medium, Rq 
the initial radius evaluated at t = to, and Ro the initial velocity, see [12, 13]. The law of motion is 



R = Ro ( l + 4^(i-t ) 
Ro 



(5) 



and the velocity 



l + 4j|(t-to) 



From equation (5) we can extract Ro and insert it in equation (6) 



(6) 



R = 



1 R A -Rt ( , R 4 - Ro^ 4 



i(t -to) Rl 
The astrophysical radius in pc as a function of time is 



1 + 



Ru 



(-) 



R(t) = ^Ro 3 (4.08 10- 6 v t (t 1 -t a ) + R )pc 



(8) 
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where ti and to are times in years, Ro is the radius in pc when ti = to and vi is the velocity in his™ units when 
ti = to. The momentum solution in the presence of a constant density medium scales as £ ' 25 . More details can 
be found in [14]. 

2.3. Classical medium with power law 

We assume that around the SNR the density of the interstellar medium (ISM) has the following two piecewise 
dependencies 

{po if R < Ro 

(9) 
po(^) d HR>Ro ■ 

In this framework, the density decreases as an inverse power law with an exponent d which can be fixed from the 
observed temporal evolution of the radius. The mass swept, Mo, in the interval < r < Ro is 

Mo = ^povr J Ro 3 • (10) 

The mass swept ,M, in the interval < r < R with r > Ro is 

M = -4r 3 po-K Hy) (d-3)- 1 

Momentum conservation in the thin layer approximation requires that 

Mv = M v , (12) 

where v is the velocity at t and vo is the velocity at t = to- The previous expression as a function of the radius is 

r 3 v (3-d) 

3r d R s - d - r a d ' V ' 

In this differential equation of first order in Ti, the variables can be separated and an integration term-by-term 
gives the following nonlinear equation Tnl 

T nl = (4i? 3 d - Ro 3 d 2 ) R - 3Ro d R 4 ~ d + R 4 d 2 
+12R 3 v t + 3i? 4 - 47? 4 rf 



+7Ro 3 vodto + Ro 3 vod 2 t 
-7Ro 3 v dt - 12R 3 vot„ - R 3 v d 2 to = . (14) 
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An approximate solution of J-nl{t) can be obtained assuming that 
3R d R 4 ~ d > -{mid - R A d 2 )R 

R(t) = 

(R 4 - d -id J R 4 ~ d (4-d) 
+ ^(4-d)v R ' i ~ d (3-d)(t-t ))^ . (15) 

Up to now, the physical units have not been specified: pc for length and yr for time are the hybrid units usually 



k m 



adopted for SNRs. With these units, the initial velocity vq is expressed in ^ and should be converted into 



This means that vq = 1.02 10 «o,i where vo,i is the initial velocity expressed in — ; in this transformation we 
have used lyr = 3.15 x 10 7 sec and 1 pc = 3.086 x 10 18 cm. The astrophysical version of the above equation in pc 
is 

R(t) = 

(i? ° 4 d " l dR ° 4 ~ d ( 4 - d ) + 3 - 402 10 ~ 7 x 

x(4-d)v 1 R s - d (3-d)(t 1 -t ))^3pc , (16) 

where ti and to are times in years, Ro is the radius in pc at ti = to and vi is the velocity at ti = to in ^j L . The 
approximate solution (15) has the following limit as t — > oo 

R(t) = C th t^z , (17) 



where 



On imposing 



we obtain 



3 Hn I 



d=^l . (20) 
a 

where a is an observable parameter defined in Section 2.1. This means that the unknown parameter d can be 
deduced from the observed parameter a. The numerical result is reported in Fig. 2. More details can be found 
in [15]. 
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Radius of SN1993J in the power law model 



o 
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Figure 2. Theoretical radius in pc as function of the time in yr as obtained by the solution of the non linear equation (14) (full 
line) which models the medium with variable density. The parameters of the simulation arc d — 2.24 , Rq — 0.031 pc, 
to — 0.049 yr and ^o,i —50 000. The astronomical data of SN 1993 J at 6 cm arc represented through the error bars 
and arc extracted from Tabic 1 in [11] 

2.4. Classical porosity 

The thin layer approximation with porosity p in classical physics assumes that only a fraction of the total mass 
enclosed in the volume of the expansion accumulates in a thin shell just after the shock front. The global mass 
between and Ro is ^irpRo where p is the density of the ambient medium. The mass that resides beyond the 
expansion is 

M = (^pRl)p . (21) 

The mass swept between and R is 

M=(^pR 3 )* . (22) 

The parameter p defined in the interval [0, 1] has been introduced because the observed radius time relationship 
in SN 1993J can be approximated by a power law dependence of the type R oc t ' 82 [11]. The theory suggests 
R oc i? 2 / 5 for the Sedov solution, see equation (1) or R oc R 1 / 4 for the thin layer approximation, see equation (5); 
therefore these two models does not match with the observations. We will see later on how is possible deduce p 
from the observed data. The conservation of radial momentum requires that, after the initial radius Ro, 

Mv = M v , (23) 
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where R and v are the radius and velocity of the advancing shock. In classical physics, the velocity as a function 
of radius is 



and introducing j3o = ^ and /3 = — we obtain 



« = «o(^)i , (24) 



P = M^f» ■ (25) 



The law of motion is 

fl(t)=(ilo 1+3F " 1 + (3 + p)«o-Rb 8, '~ 1 (*-*o)P" 1 )' fe . ( 26 ) 
where t is the time and to is the initial time. In classical physics, the velocity as a function of time is 

N 

v = „ (27) 
where 

p 

N= (iJ 1+3p_1 + (3 + p) v Ro 3p ~ (t-t )p- 1 ) 3+p v Ro 3p ~\ 
D = R 1+3p ~ 1 p + 3vo Ro 3p ~ 1 t+ vo Ro 3p ~^tp 

3 -1 3 _ 1 

-3 vo Ro p t - vo Ro p top- 

(28) 

Equation (26) can also be solved with a similar solution of type R = K(t — to) a , k being a constant, and the 
classical result is 

m= ((*+P)*Bo*''(t-t )\** (29) 



The similar solution for the velocity is 

v(t) = (3 + p)~^ p^v ^ Ro^ (t- t )~^ . (30) 

The similar formula (29) for the radius can be compared with the observed radius-time relationship reported as 

R(t) = W"*' , (31) 

where the two parameters r t s and a t s are found from the numerical analysis of the observational data, see 
Section 2.1. In this case the velocity is 

v(t) =r obs a obs t ( - a ° b °- 1) . (32) 
The comparison between theory and observation allows of deducing p: 

p = — . (33) 

1 — 0L o bs 

The theoretical solution as given by Equation (26) can be found through the Levenberg-Marquardt method 
(subroutine MRQMIN in [16]) and Fig. 3 reports a numerical example. 
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Radius of SN1993J with porosity 




Figure 3. Theoretical radius in pc as function of the time in yr as obtained by the solution of Equation (26) (full line) which 
models the medium with constant density and porosity. The parameters of the simulation are p— 18.8 , Rq — 5 10 
pc, t = 4.98 10~ 6 yr and u ,i =30 000. The astronomical data of SN 1993J at 6 cm are represented through the error 
bars and arc extracted from Table 1 in [11] 



2.5. Classical porosity and photon losses 

The emission of photons produces losses in the momentum carried by the advancing shell. This effect can be 
parametrized by the following conservation law 



(4/3 pi: R 



3\P 



v + C,(4:/3pTvY 



(R 2 - Ro 2 ) v = M vo 



(34) 



where C; (4/3p7r) p (_R 2 — Ro 2 ) takes account of the losses and C; is a constant. At the moment of writing is 
not possible to deduce the parameter Ci which controls the quantity of momentum carried away by the photons 
due to the presence of another parameter p which controls the quantity of swept mass during the expansion. The 
expression for the velocity in presence of porosity and momentum losses by a radiative process is 

, 3 

3o 

(35) 



(Rp- 1 ) 3 + CiR 2 - CiRo 2 

Fig. 4 reports how a given value of Ci modifies the velocity. The non linear equation for the trajectory in presence 
of porosity and losses is 

. 3 



Ro p 



3 + p 



^ +2 ORl +2/3 C^p + 



R p 



3+p 



3 + p 



+ 



C,R 3 C,R 2 Rp _ CiR 2 R 



3 + p 



3 + p 



3 + p 



3+p 



3+p 



- (i?o p ') /3 c(t-W) = 



(36) 
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porosity with losses 




0.02 0.04 



Figure 4. Theoretical classical velocity with porosity as given by Equation (24) (full line), velocity as given by the power law 
dependence as in Equation (32) (presence of porosity) (dash-dash line) and velocity with porosity and photon losses as 
given by equation (35) when Ci— 50 (dash-dot-dash line). The parameters of the simulation are p— 18.8 , Ro — 5 10 — 
pc, t = 4.98 irr 6 yr and w 0jl =30 000. 



An example of the modification of the trajectory as given by the parameter C; is reported in Fig. 5. 



3. Relativistic thin layer approximation 



The Newton's second law in special relativity is: 



_, dP d , . d . mnv 
F = — = — (mv) = —(- 
dt dt y ' dt 



(37) 



where F is the force, P is the relativistic momentum, m is the relativistic mass, mo is the mass at rest and v 
is the velocity, see equation (7.16) in [17]. This Section reviews the relativistic thin layer approximation in a 
medium with density regulated by a power law. A new relativistic equation of motion is introduced which works 
in a medium with constant density but adjustable swept mass. 

3.1. Relativistic power law 

In the case of the relativistic expansion of a shell in which all the swept material resides at two different points, 
denoted by radius R and radius _Ro, the previous equation (37) gives: 



M- 



M 



go 



(38) 
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Figure 5. Theoretical radius in pc as function of the time in yr as obtained by the solution of Equation (26) with porosity (full 
line) and radius with with porosity and losses as given by Equation (36) (dashed line) when Ci— 50. The parameters of 
the simulation are p=18.8 , i?o = 5 10~ 5 pc, t = 4.98 10 -6 yr and u ,i =30 000 . 



where /3o=Vo/c, /3=v/c, M is the rest mass swept between and R and M u is the rest mass swept between 
and Ro- This formula is invariant under Lorentz transformations and the initial velocity, vo, cannot be greater 
than the velocity of light. Assuming a spatial dependence of the ISM as given by formula (9), the relativistic 
conservation of momentum gives 



3 (-3 + d) sJY^W 
According to the previous equation, /3 is 



3^1 -/V 



(39) 



P = 



(-3 + flQ Porp 3 



(40) 



with Dp = 9i?o 6 /3o 2 -6Ro 6 Po 2 d 

i n 6-2 d D 2 d 

+9r R 



-6r 3 - d R a d+3 d + R 6 d 2 



n a 2 6-2 dn 2d , a a 2 3-d D d+3 

-9/30 r Ro +opo r Ro c 



In this differential equation of first order in r, the variables can be separated and the integration can be expressed 
as 

(41) 



r R 

/ ^D^dR = c (3 - d) PoRo 3 (t - to) 
Jr 



L. Zaniiietti 



Radius of SN1993J with relativistic power law 




Figure 6. Theoretical radius as obtained by the solution of the relativistic equation (41) which models the medium with variable 
density (full line). The parameters of the simulation are d — 2.52 , Rq — 0.0068 pc, — 0.14 yr ; and /3o — 0.30. The 
astronomical data of SN 1993J arc represented with vertical error bars. 



The integral of the previous equation can be performed analytically only in the cases d = 0, d — 1 and d = 2, but 
we do not report the result because we are interested in a variable value of d. The integral can be easily evaluated 
from a theoretical point of view using the subroutine QROMB from [16]. The numerical result is reported in Fig. 
6. More details can be found in [15, 18]. 

3.2. Relativistic porosity 

According to equation (37) the Newton's law in special relativity is 



F _ dp _ d_ d M r v _ d M r v 

dt ~ dt ( ' ~ dt ( yi - ^ ~ dt \/W^ ' 



(42) 



where F is the force, p is the relativistic momentum, M is the relativistic mass, M r is the mass at rest, v is the 
velocity, /3 = —, and c is the speed of light. In the case of the relativistic expansion of a shell in which the swept 
material is part of the mass contained in the advancing layer, the equation of motion is 



(43) 



where Po—Vo/c, Ro is the initial radius, and R is the temporary radius. The velocity of a relativistic expanding 
shell is 

PoRo 3/p 



^R^' p Po 2 + R 6 /p - R e /pPo 2 



(44) 
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Figure 7. Theoretical relativistic velocity for SN 1993J as given by Equation (44) (rclativistic porosity-dashed line), theoretical 
classical velocity as given by Equation (24) (classical porosity -full line), velocity as given by the power law dependence 
as in Equation (32) (dot-dash-dot-dash line). 



Fig. 7 shows the classical and relativistic behaviors of the velocity as a function of the radius R. In this differential 
equation of the first order in r the variables can be separated and the integration can be expressed as 



[ R Jr 6/p Po 2 + i? 6 / p - R 6 /Pj3 2 dr = j3 R 3/p (t - t )dt 



(45) 



The integral of the previous equation can be expressed analytically but the result is complicated. The integral 
can be easily evaluated numerically using the subroutine QROMB from [16]. The numerical result is reported in 
Fig. 8. An interesting property of the Gamma-ray burst (GRB) is that the frequency of observation at which it 
starts to be visible decreases with time. An example of such behavior for GRB 050904 in the BAT, XRT, J, and 
I bands, can be found in Table 1 of [19] . The numerical analysis of the data as reported in Fig. 9 gives 



v = 1.03 10 25 t(sy 2 - 27 Hz 



(46) 



where t has been chosen as the value of time at which the flux at the chosen frequency is maximum. The interval 
in time at which the GRB becomes visible in different astronomical bands is here correlated with the velocity of 
the shell's expansion, see Fig. 10. 

4. Conclusions 



We have deduced two new laws of motion for an advancing shell assuming that only a fraction of the mass which 
resides in the surrounding medium is accumulated in the advancing layer. According to the present laws of 
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Radius of SN1993J with relativistic porosity 




Fi gur6 8. Theoretical radius in pc as function of the time in yr as obtained by the solution of the relativistic Equation (45) 
(relativistic porosity) (full line). The parameters of the simulation arc p— 18.3, Rq — 6 10 -5 pc, to — 5 10~ 8 yr, and 
{3 =0.1. The astronomical data of SN 1993J arc represented through error bars and are extracted from Table 1 in [11] 




Log (t [sec]) 



Figure 9. Bands from the radio through the IR/optical to X-ray and BAT at which the expansion becomes optically thin. The 
energy bands are expressed in Hz and the y-axis as the logarithm of the frequency. The x-axis reports the logarithm of 
the time in s. The data are extracted by the author from Table 1 in [19]. 
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velocity of GRB with relativistic porosity 
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Figure 10. Theoretical relativistic velocity for SN 1993 J as given by Equation (44) as function of Log(timc) in seconds. The 
parameters of the simulation are p— 18.3, Rq — 2.25 10~ 8 pc, to — 5 10 — 8 yr, and /3o —0.33. 

mechanics the analysis has been split into classical and relativistic cases, see Equations (26) and (45). In the 
classical case a similar solution has been found, see Equation (29). The momentum carried away by the photons 
can be modeled in the classical case by introducing a reasonable form for the photon's losses , see ( 36). 
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